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Data  taken  from  the  second  generation  of  the  Naval 
Postgraduate  School’s  ocean  floor  aid  land  geomagnetic  data 
collection  system  is  analysed  to  determine  the  wave  nature 
of  Geomagnetic  Hicrcpulsat ions  in  the  .35  to  13  Hz  range. 
The  measurements  were  made  in  Honterey,  California  during 
June  to  September  1982.  Three  of  these  experiments  were 

twenty- four  hours  in  length.  During  two  such  exoerim*nt- 
land  and  sea  data  were  transmitted  to  a  central  point  where 
the  signals  were  recorded  simultaneously.  Geomagnetic  fluc¬ 
tuations  were  analyzed  to  determine  various  factors  charac¬ 
terizing  the  waves.  These  factors  include  the  power  spectral 
densities,  coherences,  Stores  parameters,  orientation  angle, 
degree  of  polarization  and  ellioicity.  The  software  devel¬ 
oped  in  this  research  provides  an  ettensive  capability  for 
lysis  of  geomagnetic  noise,  ions  term  data  collection  can 
utilize  it  to  establish  trends  in  octh  the  spatial  and 
temporal  properties  of  naturally  oc c n ring  EL?  electromag¬ 
netic  fields. 
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I.  [  NTROgtJCTIDS 


This  thesis  is  part  of  the  long  term  effort  of  the  Naval 
Postgraduate  School  to  collect  and  interpret  EL? 
elctromacnetic  noise  on  tie  sea  floor  and  from  a  nearby  land 
station.  The  project's  oojectiva  is  to  study  and  interpret 
signals  in  the  0.05  to  13  dz  freuuencv  ran a®.  Z t  “mohasize? 


not  se 


rm  cert  a  nee  or  octant  mg  measurements  r : 

:  on  the  sea  fleer  and  or.  land  over  a  eerie 


:eoaa  jr.etio 


several 


years. 

The  specific  obectives  of  this  thesis  are  twofold;  to 
calibrate,  install,  and  operate  tie  land  data  collection 
system;  and  to  develops  a  computer  program  to  analyze  the 
data  recorded  to  provide  easily  utilized  interpretive 


A.  IT.tZ  V  A  r.  I  AT  I  CMS  07  7  72  GEOMAGNETIC  ?~ZLD 

The  earth's  magnetic  field  has  oeen  used  for  centuries 
for  navigation  at  sea  aid  for  survey  direction  on  land. 
These  common  uses  are  dependent  on  the  static  magnetic 
field.  The  static  field  can  be  generally  characteriz ed  as  a 
magnetic  dipole  at  a  11.5  degree  inclination  to  the  axis  of 
rotation  of  the  earth.  It  shows  only  minor  diurnal  varia¬ 
tions  ir.  direction.  Its  magnitude  varies  greatly  over  the 
earths  surface  but  usually  is  in  the  range  of  43,330  to 
70,  000  nano-Ieslas  (nT).  It  Monterey,  California  the  static 
field  has  a  valu<=>  of  48,633  nT. 

It  has  only  been  since  the  discovery  of  magnetic  induc¬ 
tion  that  the  underlying  variations  In  the  static  field  have 
become  well  known.  These  variations  are  due  to  many 
sources.  The  sources  include,  but  are  not  limited  ko.  solar 
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magnetic  storms. 


activities,  lunar  daily  variations, 
derstorms  in  the  earth's  atmosphere, 
detonation  of  nuclear  davit  as.  One  o: 
to  this  is  the  interaction  of  the  cha 
sua  with  the  earth's  magnetosphere.  r 
in  duration,  amplitude,  aid  the  frega 
within  the  frequency  range  of  this 
variations  are  characterized  as  micro 
Gecmaar.et  ic  micr  coul  si  t  ior.s  are 


earthquakes. 


t  h  u  r.  - 
and  the 


the  major  contributors 
rged  particles  from  the 


hase  events  vary 

widely 

ency  range  they  e 

f  £  o  r*  — 

system  the  most 

common 

pulsations. 

defined  tc  os 

■ic*  ua- 

tie  ns  in  th«  earth's  magnetic  field 


0.2 

seconds  to 

1  0  mi 

r.ute  s.  They 

vary 

r:T 

tc  as  high 

as  a 

few  tens  of 

nT. 

tha 

interaction 

of 

groups  of  c 

harged 

ir.  the  earth's  magnetosphere  as  hydro 
waves  inturn  create  the  soiree  curre 
icasphere),  for  the  micropulsation  ? 
are  divided  into  two  sub categories: 
nuous.  The  abr aviations  for  these 
continuous  and  pi  for  pilsaticn  L 
difference  between  the  two  is  that 
have  a  semi- sin  usoiaal  nature  while 
show  irregularities  in  both  amplitud 
period  ranges  and  average  amplitude  i 
listed  in  Table  I,  [  Hef .  1]. 

The  overwhelming  majority  of  t 
frequency  range  are  due  to  solar  a 
"Schumann"  resonances,  will  also  be  a 
of  6-8  Hz.  These  resonances  are  due  t 
in  the  concentric  spherical  cavity 
surface  and  the  lower  region  of  the  i 


which  have  per:  * -  from 
in  amplitude  f:  .001 

Tney  are  attributed  tc 
particles  propagating 
n agnatic  waves.  These 
ats,  (primarily  ir.  the 
ield.  Micropulsations 
irregular  an d  cor. ti¬ 
ers  pc  for  pulsation 
rregular.  rue  primary 
tie  pc  micropulsations 
tie  pi  micropulsations 
e  and  frequency.  The 
or  th»se  pulsations  are 

he  pulsations  ir  this 
ctivity.  In  addition, 
oparer.t  in  the  vicinity 
o  lightening  transients 
between  the  earth's 
onosphere,  "Bef.  2]. 


The  sensor  system  currently  employed  consists  of  a  thrsa 
coil  land  system  deployed  at  La  Mesa  Village  nousing  area, 
approximately  2  km  from  tie  Naval  Postgraduate  School  data 
collection  center,  and  a  two  coil  sea  system  deployed  or.  the 
bottom  cf  Monterey  Bay.  Its  distance  to  the  Naval 
Postgraduate  School  (NPS)  varied  from  2  to  4  km.  The 
original  intention  of  the  project  <as  to  deploy  the  land 
system  approima tel y  22  miles  from  tne  Naval  Postgraduate 
Scicci  at  Chew's  Ridge.  This  prove?  impossible  due  to  data 
transmission  difficulties.  See  Apoer.dix  3  for  additional 
information  on  this  site. 

Each  sensing  coil  is  a  continuously  wound  coil  antenna 
made  up  of  approximately  54  60  turns  of  18  gauge  copper  wire. 
The  average  sensing  area  of  each  coil  is  3.0825  square 
meters  and  its  dimensions  are  as  depicted  in  Figure  1.1. 
The  sea  system  has  its  coils  mounted  orthogonally,  usually 
in  the  horizontal  plans.  The  land  system  coils  are  affixed 
to  a  wcccer.  mount  that  holds  them  in  ar.  orthogonal  configu- 


sea  system  cannot  be  determined  when  deployed.  The  land 
system  is  oriented  witn  tne  X-coil  in  the  magnetic  ncrtn- 
south  direction  ,  ?-coil  In  the  magnetic  east-west  direc~ion 
and  the  Z-coil  in  the  vertical. 

The  two  systems  were  designed  to  be  as  similar  as 
possible.  The  coils,  preamplifier,  signal  conditioner  and 
pulse  code  modulator  (PCM)  are  identical  in  design  and  were 
all  calibrated  to  insure  tne  production  of  valid  data.  Th= 
technical  description  of  iLl  of  these  items  may  be  found  in 
[Ref.  3]  .  The  land  system  is  deployed  as  shown  in  Figure 
1.3.  The  pre- amplifiers  are  battery  powered  and  located 
within  one  foot  of  the  colls.  The  coil  mount  and  preampli¬ 
fiers  were  positioned  approximately  50  meters  from  tne 
signal  conditioners,  and  connected  by  three  coaxial  cables. 


c 


Figure  1.2  Sensor  Noaatmg  Block 

The  signal  coni i tioners ,  PCM  boar!  and  transmitter  are  al 
housed  in  a  small  wooden  onilding  and  powered  by  D.3.  pcwa 
supplies.  This  building  is  the  closest  source  of  alternating 
current. 

After  the  data  has  been  encoded  oy  the  PCM,  it  is  tran¬ 
smitted  by  7 HF  radio  data  link  to  the  Naval  Postgraduate 
School  for  recording. 


C.  DATA  COLLECT  ION  AND  PROCESSING 

Land  and  sea  data  are  simultaneously  received  an 
recorded  on  separate  channels  of  an  analog  tape  recorder 
Each  tape  holds  60  to  120  ainutes  of  data.  These  tapes  mus 


Figure  1.3  Land  Data  Collection  System 

be  decoded  and  re-recorded  onto  digital  tape  suitaole  far 
the  main-frame  computer  ac  the  Chum  Computer  Center.  Each 
digital  tape  contains  93  ninutes  of  Land  cr  sea  data.  Each 
of  these  digital  tapes  is  then  processed  through  the 
ccnputer  program,  (See  Ciaptar  3,  Section  C.l ,  to  provide 
the  piers  of  the  various  interpretive  factors.  See  3eferer.ee 
3  for  the  details  of  the  recording,  decoding  and  re-re¬ 
cording  on  diaital  tape. 


II.  SYSTEM  CALIBRATION 


A.  INTRODUCTION 

The  triaxial  coil  magnetometer  system  is  diagramed  in 
Figure  2.1.  The  lata  sensed  oy  this  magnetometer  system  is 
giver,  by  Vx  (t)  ,  Vy  ( t)  and  1£  (t»  ;  -.he  cutout  voltages  gener¬ 
ated  by  the  input  magnetic  fields.  Tie  sys-.em  transfer  func¬ 
tion  must  be  determined  sc  that  the  magnetic  field 
components  can  be  removal  fr on  the  oatpu-  voltage.  In  the 


tine  domain  the  relationship  between 

the  output  vol 

tage  and 

the  magnetic  field  is  a  convolution 

lot  sgral. 

v;  (M  =  /h ;  (t-t')o;  <t«.  if 

(1> 

where  i=x,y,z.  It  is  more  convenient 

to  express  Sgn. 

(i>  in 

the  frequency  domain. 

Vj  (<*»)  =H;  (u<)  3  j  (ut)  .  (2) 

The  transfer  function  i («-u)  can  easily  be  determined  in 
terms  cf  the  ratio  of  output  to  an  input  test  signal  cr  by 
analysis  of  the  measures  an  t  system's  properties. 

The  next  section  discusses  the  components  cf  the  system 
and  attempts  to  develop  i  system  nodal.  The  last  section 
describes  the  calibration  measurements  and  the  experimental 
determination  of  the  transfer  function. 

B.  SYSTEM  DESCRIPTION 

The  individual  transfer  function  for  the  sensing  coil 
is: 

VcM  =  j^KB  {*/)  (3) 


U 


I 


I 


Figure  2.1  Systei  Coaponent  Configuration 

where  Vc  (»)  is  the  spectrum  of  tae  voltage  directly  out  of 

the  coil  and  K  is  a  gain  factor. 

The  r.ex-1’  component  in  the  system  is  the  preamplifier  low 
pass  filter.  The  preampLLfier  introduces  a  gain  of  1000  so 
the  voltage  signal  can  be  sent  to  remote  instrumentation. 
The  low  pass  filter  is  a  five  poLe  Chebyshev  with  a  3dB 
cutoff  at  20  Hz.  The  filter  is  primarily  used  to  remove  50 
Hz  signals  from  the  system  and  to  prefilter  the  data  before 
if  is  digitized.  An  additional  pole  is  introduced  because  of 
a  capacitive  coupling  between  the  second  and  third  stages  of 

the  filter.  Figure  2.2  shows  a  pole-zero  diagram  of  the 

preamplifier.  The  transfer  function  tafces  the  form: 

H(j«0  *K/(j«+a)  (jw+b)  (j«w*o-  jd)  (jw  +  c«-jdi  (j^+e-jfi  (ju/+  =  +  jf)  (t*i 
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where  K  is  the  gain  of  tie  preamplifier  and  the  coeficisnts 
a  -  f  are  determined  from  the  pcle-zeco  diagram. 

The  preconditioner  aas  additional  gair.  of  approximately 
1030  and  has  zener  diodas  to  limit  tie  input  voltage  to  the 
PCS  board.  Typically  tha  cutoff  is  around  7  to  3  volts. 

The  magnitude  and  phase  responss  of  the  preamplifier- 
filter  and  preconditions r  are  shown  in  Figures  2.3  through 
2.3.  The  input  test  signal  generating  these  curves  had  flat 
magnitude  and  phase  ciaraoter  as  sh-'*"  *  *«  v-Ti-as  ~>  Q  = 

figures  then  represent  tns  actual  t 
pra amplifier- filter  pra conditioner 
channels  of  land  and  ocean  sensor  sy 


ar  fun 


or.s  of  the 


s  show 


1  -hi 

t  the  channels 

are  welL  matched  u 

magnit 

'  Tha 

phase  response 

is  linear  and  vir 

tualiy 

out  ou  . 


every  case.  This  is  an  important  ?o 
of  the  transfer  function  is  not  removed 
voltage  signal.  However,  since  tha  data  analysis  is  in  terms 
of  cover  spectral  densities  and  coherences,  this  fact  has  no 
oearing  cn  the  results. 

The  total  transfer  function  of  tna  system  is  tha  product 
of  Egn.  (3)  and  (4)  .  Based  on  thesa  formulas  and  tha  figures 
tha  resulting  transfer  function  is  predicted  to  be  fairly 
linear  at  the  low  frequency  end. 


C.  DETERMINATION  OE  THE  5?STSM  TRANSFER  FUNCTION 

The  magnitude  of  the  system  transfer  function  was  deter¬ 
mined  experimentally  through  ta?  uss  of  a  Helaaoltz  coil,  a 
source  of  sinusoidal  currant  and  a  spectrum  analyzer.  The 
experimental  set  up  is  illustrated  in  Figure  2.10.  A  ona 
turn  Helmholtz  coil  with  a  0.51  meter  radius  and  a  0.51 
meter  seperation  between  coils  was  isai  tc  provide  a  cali¬ 
brated  magnetic  field  to  tne  sansi.ng  coil.  The  sensing  coil 


IB 


Pigura  2.6  Response  af  the  Dcean  1  systen 


was  placed  concentric  tc  the  Helmholtz  coil,  and  cantered 
vertically,  to  ensure  a  uniform  magnetic  field  about  the 
sensing  coil.  The  magnetic  field  in  this  region  of  the 
Helmholtz  coil  is  directed  along  its  axis  and  is  given  by 

[  Ref.  4  1  : 

3  =>4,8  N 1/5  y^B.  (5) 

to  terms  on  the  order  of  (L/R),  where  L  is  the  linear  dimen¬ 
sion  of  the  region,  N  is  the  number  of  turns,  ?.  is  the 


i.  Tr.e  1  riving  currer. 
erazcr,  Nodel  Number 


by  a  tfavetek  signal 
ent  was  measured  bv 


ding  the  voir  age  across  a  996  ohm  resistor  in  series  with 
Helmholtz  coil.  A  spertrum  analyzer  (HP-3582A)  was  used 
measure  the  system  outpu*  tc  avoid  a::y  spurious  signals 
ated  by  harmonics.  By  repealing  the  measurements  at  a 
number  cf  different  frequencies  and  drive  currents,  the 
svstem  transfer  function  could  be  measured. 


The  actual  calibration  was  conducted 
villaae  site  usir.u  •‘■he  setuc  illustrated 


i  u  r  s 


2.  19. 


The  Helmholtz  ceil,  sensing  ceil  and  preamplifier  were 
located  approximately  40  meters  from  the  shack  which  housed 
the  remainder  of  the  experiment.  4ith  the  Helmholtz  coil 
disconnected  the  background  magnetic  noise  of  this  location 
was  repeatedly  examined  to  ensure  accurate  readings.  Figure 
2.11  shows  a  typical  background  noise  plot  from  the  spectrum 
analyzer  (the  transfer  function  has  not  been  removed).  The 
distinct  features  are  Schumann  resonances,  otherwise  the 
background  was  low  enough  and  stable  enough  to  adequately 
calibrate  the  system.  Only  at  the  Schumann  frequencies  was 
the  background  subtracted  from  the  neasured  output  voltage. 
Calibration  runs  were  conducted  from  3.95  Hz  to  20.9  Hz  with 
applied  calibrated  magnetic  field  of  0.02,  0.2  and  1.0 


Figure  2.10  System  Calibration  Experiment 


nancteslas.  Figures  2.12  through  2.16  show  the  resulting 
calibration  curves  of  the  land  and  ocean  magnetometer 
systems.  The  true  transfer  function  is  represented  by  the 
0.02  nT  field.  The  larger  applied  fields  show  saturation  of 
the  system  at  the  protective  zener  iiodes  in  the  precondi¬ 
tioner.  At  the  low  frequency  end  all  the  systems  are  closely 
matched.  At  1  Hz  an  output  voltage  of  1  v  indicates  a 
measured  magnetic  field  of  0.4  nT.  A  1  nT  field  is  saturated 
above  2  Hz. 

The  comDutation  of  the  equations,  representing  the 
magnitude  of  the  transfer  functions,  requires  fitting  the 
data  to  the  theoretical  transfer  fanctior.  of  the  previous 
section.  This  would  be  desirable  if  the  phase  character  of 


the  sys-em  transfer  function  need  be  estimated.  However, 
since  the  primary  area  of  interest  is  in  the  lower  frequency 
region,  and  since  the  phase  response  is  net  needed  if  the 
various  channels  maybe  assimed  identical,  it  was  decided  to 
model  the  transfer  functions  in  terms  of  linear  segments. 
The  equations  derived  in  this  manner  have  a  maximum  error  of 
If  in  the  range  of  0  to  1 0  Hz  and  5s?  in  the  10  to  20  Hz 
range  for  the  magnitude  of  the  transfer  funtior.. 

In  the  actual  computation  each  curve  was  broken  into 
convenient  segments  that  could  be  represented  by  a  straight 


line.  The  range 

of  the  segm 

er.ts  were 

0  tc  5  Hz, 

5  tc  7.5  H 

7.5  to  10  Hz, 

10  to  15  iZ 

and  15  to 

20  Hz. 

Ih e  str ai g! 

line  parameters 

were  compi 

ted  for  ? 

ach  segment 

to  yield 

series  cf  five 

equations 

tc  model 

each  ser.si 

n  g  coil  a 

amplification  system.  See  Table  II  for  the  individual 
transfer  functions. 
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Figure  2.11  Background  3ao«agnetic  Moi s?  f  (Hz)  05/20/82 
15  13  hrs. 
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Figure  2.14  Laai  Systea  Z 
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Pigure  2.15  Ocean  Systea  X-3ail  Calibra 
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•  2 ill  2V ALO A I T D N 

A.  POWER  SPECTRAL  DENSITY  AND  ALIASING 

The  first  step  in  the  analysis  of  the  data  collected  -was 
the  determination  of  the  pc  we  r  spectral  densities  (PSD's) 
for  each  channel  of  the  land  and  sea  lata.  This  was  accom¬ 
plished  * hr c ugh  the  use  of  the  computer  program  developed  by 
J.  a.  Schweiger  fRef.  5],  After  minor  modification  of  this 
program,  approximately  twenty  hours  of  real  time  data  was 
evaluated.  Samples  of  the  PSD's  are  presented  in  Figures  3.1 
to  3.5. 

The  overall  value  of  the  PSD's  is  ir.  agreement  with 
previous  work  conducted  in  this  area.  The  PSD's  have  a 


decline  cf 

appr 

cximatelv  -6 

d3  per  octave 

from  0  to 

5  Hz 

where  they 

t  er.  d 

to  decrease 

their  rate  of 

decline. 

Night 

PSD's  tended 

+  £ 

be  3  to  7  13 

below  tie  dayt 

This 

coincides  with  the  findings  of  NcDe vi * t  and  Homar.  (1930),  :■]. 
W.  Beard  (1981)  and  A.  C.  Fraser- S mi th  and  J.  L.  Buxton 

(1975),  [Pef.  6]. 

One  feature  of  the  PSD's  that  differed  from  previous 
PSD's  is  the  spike  at  4  HZ  ir.  the  land  data.  This  was  found 
to  be  "aliasing”  of  the  53  Hz  power  grid  component  by  the 
PC!!  sampler.  This  experiment  was  conducted  with  a  sample 
rate  of  22  Hz;  which,  by  the  sampling  theorem,  should  give  a 
clear  view  of  the  frequency  range  of  0  to  10  Hz,  when 
coupled  with  the  20  Hz  cut-off  filter  in  the  amplifier, [  Ref. 
7%  The  60  Hz  component  cf  the  geomagnetic  field  was 
measured  by  F.  W.  Clayton  (  1973)  end  found  to  be  4  (r.T)  /Hz, 
which  is  10  higher  than  the  surrounding  data  ‘Ref.  8].  The 
20  Hz  filter,  built  into  the  amplification  system,  would 
normally  eliminate  frequencies  in  this  range  but  because  of 
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the  great  differences  in  magnitude  this  does  net  happen. 
When  the  32  to  64  Hz  range  is  convolved  onto  the  0  to  32  Hz 
range  the  60  Hz  component  will  appear  at  4  Hz.  This  destroys 
the  usefulness  of  the  data  in  a  snail  region  about  this 
poi nt . 

Another  point  of  interest  is  illustrated  by  the  vertical 
component  of  the  sea  data  recorded  on  September  16,  1982, 
Figure  3.5.  During  this  experiment  the  Y-coil  was  oriented 
in  a  vertical  rather  than  the  more  normal  horizontal  direc¬ 
tion.  Examination  of  the  da~a  shows  11  to  15  13  decrease  in 
the  PSD.  This  is  in  agreement  with  the  theory  that,  oecause 
of  the  ccean-a ir  interface,  geomagnetic  waves  propagate 
vertically  in  the  sea,  this  eliminating  this  component. 


S.  STOKES  PARAMETERS  AMD  THEIR  INTERPRETATION 


A  ft  sr 

the 

collection  of 

the 

data  during 

June  to 

September 

of  198 

2,  the  problem 

of  how 

to  interpret 

it  became 

a  orime 

concern 

.  D urine  the  1 

after 

o a r t  of  this 

c  —  r  1. 0  ^  2. 

computer  program  was  developed  to  determine  the  Power 
Spectral  Densities  (?SD's>  by  J.M.  Sohweiaer,  October  1982, 


[Re-- 


5]. 


i 


This  program  provided  the  basic  structure  tc  read  the 
data  points  off  of  magnetic  data  tape  and,  through  the  use 
of  a  Fast  Fourier  Transform  (FFT)  routine  and  transfer  func¬ 
tions,  to  compute  the  geomagnetic  fluctuations  in  the 
frequency  domain.  This  program  provided  the  basic  format 
which  was  expanded  and  modified  to  produce  other  parameters. 
The  derivations  of  the  equations  for  these  parameters  and 
their  interpretation  comprise  the  remainder  of  this  section. 

The  first  parameter  to  be  discussed  is  the  coherence 
between  the  different  channels  on  land  and  sea.  Coherence 
gives  a  measure  of  the  "singleness''  of  J-he  source  at  the 

"Singleness"  in  this  context  refers  to 


various  rreguer.cies 


PSD,  Vertical  CDipon 


whit  percentage  cf  the  intensity  at  a  frequency  can  be  saii 
to  come  from  a  single  source.  To  compute  the  coherence  one 
first  calculates  the  cross  power  spectral  densities  between 
the  channels.  For  a  given  pair  of  channels  X  and  Y ,  the 
cross  power  spectral  density  (S  )  is  defined: 

N 

Siy(fj)=1/N^X(fj)  Y(fj) 
n=  1 

In  this  equation  u  is  the  number  of  blocks  of  data  sent 
through  the  FIT  and  averaged;  X  ( fj )  and  Y(f;)  are  the  FFT’s 
along  the  appropriate  axis.  The  co  herence  (Vty\  is  defined 
as: 

S „(?j  ) 


(S*  C-J  )  Sr(f;  )  ) 

For  the  land  data  coherences  were  also  determined  for  the 
Y-3  plane  and  the  2-X  plane.  Figures  3.7  and  3.3  show 
typical  land  and  sea  coherences,  with  all  coils  oriented  in 
the  horizontal  plane.  Note  that  on  the  land  the  coherence 
stays  very  close  to  its  maximum  value  of  one  until  approxi¬ 
mately  7-8  Hz  and  then  starts  to  decrease  while  in  the  case 
cf  the  sea  data  there  is  more  structure  at  the  lower 
frequencies.  This  type  of  behavior  is  typical  of  most  of  the 
data,  although  the  exact  frequency  of  the  decline  varies. 

In  addition  to  coherence,  Stokes  parameters  may  be 
computed  from  the  auto-spectral  densities  and  cross-spectral 
densities.  They  can  be  used  to  interpret  the  geomagnetic 
waves  cn  a  bi-Linear  basis,  a  bi-Circalar  basis  or  an  ellip¬ 
tical  basis. 

Although  other  notations  are  utilized  for  the  Stokes 
parameters,  the  most  common  notation;  sC ,  si,  s2,  s3;  will 
be  employed  here  [Ref.  9].  The  general  form  of  ‘he  Stokes 
parameters  are: 
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Figure  3.7  Laad  Coherency  X  Vs  Y  Channels 


log  Hz 


sO=<alS+<a2Z> 
s1=<a1Z>-<a2*> 
s2=2<a1 a2  cos £  > 
s  3=2  <a 1  a2  sin£  > 

where  al  and  a2  are  the  instantaneous  amplitude  of  the 
orthogonal  components  of  the  magnetic  vector  and  £  is  their 
phase  difference.  In  terms  of  the  auto-  and  the  cross-spec- 
*ral  densities  the  Stokes  parameters  oar.  be  written: 

s0=2  (5*  *Sy) 
s  1=2  -Syl 
s2=4  3e(^y) 
s3=4la(%y) 

It  must  be  noted  that  the  Stokes  parameters  were  developed 
originally  to  describe  monochromatic  light  and  can  only  b= 
apclied  to  cur  problem  after  the  data  has  been  transforms* 
to  -he  frequency  domain. 

The  bi-Linsa  r  basis  interpretation  of  the  Stokes  parame¬ 
ters  can  be  seen  by  examining  a  superposition  of  two  orthog¬ 
onal  linearly  polarized  waves  [Ref.  10].  Computing 
Stokes  parameters  in  this  basis  yields: 

s0=  1/2  (E*  *E*  ) 
s  1  =  1/2  (B^-bJ) 
s2=  By  ByCO  s  (9f  -9y  ) 
s3=  Bx  By  sin  (Sj,  -9r  ) 

Another  basis  useful  for  interpretation  is  the 
bi-Circular  basis.  This  basis  is  represented  by  the  superpo¬ 
sition  of  two  counter-rotating  circularly  polarized  waves. 
In  this  basis  the  Stokes  parameters  are: 


s0=3%  BX 


The  important  point  to  none  from  these  different  basis 
secs  is  that  for  a  given  wave  form  (linear,  elliptical,  or 
cicular)  the  Stokes  parameters,  as  computed  from  the  spec¬ 
tral  densities,  are  the  same  in  any  basis  in  which  they  are 
computed.  This  can  be  easily  visualized  by  taking  a  circu¬ 
larly  cr  linearly  polarized  wave  and  determining  the  value 
of  the  Stokes  parameters  in  each  basis.  It  seems  that  these 
interpretations  are  most  easily  visualized  in  the  elliptical 
basis.  Consequently,  in  addition  to  the  Stokes  parameters 
the  angle  <*c  was  computed  to  make  the  interpretation  easier. 
This  ar.ale  was  calculated  from  the  Stokes  oarameters  si  and 


s2  by  the  equation: 


o<.  - 1/2  arctan  (s2/s1 ) 


Examples  of  the  results  of  the  calculations  of  the  Stokes 


oarameters  and  the  ancle  oc  for  land  ar.d  sea 


lata  ara 


displayed  in  Figures  3.10  to  3.19.  note  that  si,  s2,  and  s3 
are  normalized  by  dividing  by  sO.  In  addition,  for  the  land 
data,  the  same  graphs  were  produced  for  the  y-Z  plane  an 
the  Z-X  plane. 

The  use  of  Stokes  parameters  is  based  cn  aoncchr cmatic 
ar.d  polarized  waves.  To  determine  the  applicaDility  to  our 
data  of  the  Stokes  parameters  the  iegree  of  polarization 
must  be  found.  For  a  perfectly  polarized  monochromatic  wave 

s0l=s1*  «-s2l  +S32-. 


If  the  wave  is  partially  polarized  then  sO  becomes  greater 
than  the  sum  of  the  squares  of  the  other  three  [REF. 9].  To 
determine  how  close  to  the  polarized  case  the  data  lies,  the 
degree  cf  polarization  (P)  can  be  calculated.  It  is  defined 
by, 

(5  l"*  +32*  +s3a> 


M* 


Figure  3.17 
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If  P  is  one  then  the  wavs  is  perfectly  polarized  and  any 
analysis  based  on  the  Stokes  parameters  will  be  completely 
accurate. 

For  values  of  P  less  than  one  the  analysis  is  corres¬ 
pondingly  less  valid.  However,  factors  that  depend  only  on 
si,  s2,  and  s3,  such  as  the  angle  and  the  ellipticicy,  are 
net  effected  by  low  values  of  P.  Ia  such  regions  the  r.cise 
level  is  much  higher  and  the  data  is  non  easily  interpreted. 

The  land  and  sea  data  vary  in  their  values  cf  ? 

depending  cn  the  time  of  day  and  otner  unknown  factors.  In 
general  they  both  shew  a  high  degree  of  polarization  through 
3  Hz  while  above  this  frequency  a  marked  drop  off  occurs. 
If  was  also  noted  that  in  the  land  data  the  degree  cf  polar¬ 
ization  is  virtually  identical  in  all  three  planes.  This  is 
aoDarer.tly  due  +  he  fact  -.lit  the  polarization  was  not 

oriented  along  any  cf  the  three  sensing  axis  and  therefore 
the  proportion  cf  the  total  intensity  that  resulted  from 
polarized  waves  remained  relatively  constant.  (See  Figures 
3.20-  3.  23) 

Cne  final  tool  cf  interpretation  of  the  data  is  the 

ellipticity  of  the  wave  (E)  .  The  sllipticity  is  defined  by 

the  majer  and  minor  axis  of  the  ellips e , [ Ref .  3]. 

Specifically  E  is  defined  by: 

fan  3= 

This  car.  also  be  written  in  terms  oe  the  Stokes  parameters 
as: 

X  2  i  ^ 

E=1/2  aresm  (s3  /  ( s  1  +  s2  -t-s3  )  ) 

where  E  varies  fromfT/4  to  -  7^4 .  The  interpretation  of  E  is 
illuminated  by  considering  its  dependence  on  s3  and  the 
elliptical  basis.  If  E  is  zero  then  the  minor  axis  is  zero 
which  means  the  polarization  of  the  wave  is  linear.  At  E=^/u 
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log  Hz 


gure  3.22  Degree  of  Polarization,  Land,  Z-X  Plat 


log  Hz 


Pigars  3.25  Ellip*i=ity,  3h,  X-Y  Plm 


the  polarization  is  circular  ana  is  rotating  in  a  right- 
handed  sense,  and  an  E=- ^/4  the  wave  is  again  circular,  bu* 
with  a  left-handed  rotation.  In  between  these  values  cf  5, 
the  polarization  of  the  wave  is  an  ellipse  with  the  appro¬ 
priate  handedness  of  the  polarization  depending  on  whether  E 
is  positive  or  negative.  The  smaller  the  absolute  value  oc 
the  ellipticity  the  narrower  or  acre  linear  the  ellipse 


oeccmes. 


1  r.  e  s  e  w  ?.  v  -s  tor 


From  ?igur-s  3.2'u  and  3.  If,  -he  fora  of  -hese  w  ?.v~s  for 
both  lard  and  sea  data  is  generally  a  narrow  ellipse  with 
varying  left  and  right  handedness  in  the  areas  of  high 
polarization  while  appearing  to  be  nore  elliptical  out 
these  areas  [Ref.  9], 

The  final  result  of  this  investigation  is  that  thr 
the  computation  of  the  Stokes  parameters  the  nature  of 
geomagnetic  wave  can  be  determined.  Specifically,  tc 


e — ipticiry  ana  tne  orientation  anp^s  oc  . 


j.  .1  -  -  e  u  _  -J  -r 


pc^-ar  izar  ion  is  a.so  nee  a  ea  tc  es t a  c_isn  “he  oegree  of 
validity  of  the  interpretation.  Naval  Postgraduate  School 
Techr.ical  Reoort  *NPS  61-33  -005  ,  contains  additional  plots 
of  all  ■‘•he  factors  described  in  this  section. 


C.  COMPUTER  PROGRAM  DE7  EL  0  PEMENT  AND  DATA  PROBLEMS 

As  noted  earlier  the  ability  to  process  large  amounts  of 
data  in  a  timely  manner  was  first  put  to  use  in  this  project 
in  September  1932.  The  original  program  has  been  expanded 
to  compute  all  of  the  factors  listed  in  the  previous 
section.  It  consists  of  a  main  program,  a  FFT  subroutine, 
and  a  read  subroutine.  As  noted  before  the  FFI  subroutine 
transfers  the  data  into  the  frequency  domain  and  the  read 
subroutine  reads  the  PCM  data  from  tne  digital  computer  tape 
and  places  it  in  the  input  array.  All  changes  and  additions 
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male  were  in  the  main  program,  therefore  only  the  main 
program  is  included  in  Appendix  A.  This  program,  Lanaxyz,  is 
used  to  evaluate  only  tae  land  data.  The  version  of  the 
program  used  to  evaluate  the  sea 
Landxyz  with  the  following  exceptio 
functions  are  used;  because  the  s 
coils,  only  one  plane  of  interpretiv 
Although  the  basic  program  prov 


rrc_ude  me  ether  calculations  create!  severe^  proexems.  me 
Job  Control  Lap.  au  age  (JCL)  provides  the  guidance  tc  tha 
operating  system  of  the  computer  *:o  tell  it  hew  to  run  batch 
iebs.  Specifically  what  input  ar.i  output  devices  to  use  and 
hew  much  memory  to  allocate  *o  the  program  during  each  step 
of  its  execu-ion.  (In  the  program  the  JCL  can  be  identified 
by  the  // ' s  at  the  start  of  each  line  of  this  code.)  D urine 


or  *  n  —  •  n  e  it  duh  .  3  of  ?).  9  3  0  r  v 


-  -  -  u  - 


L ;;;  X ,  GO  and  PLOT  steps  in  the  execution  of  the  program  ha  v 
all  exoar.  iei.  The  current  JCL  is  caoable  of  handling  up  t: 
two  hours  of  data  at  32  Hz  or  one  near  of  data  at  6u  Hz. 

Sea  and  land  data  are  initially  recorded  on  analog  *ap- 
ar.d  then  decoded  ar.d  re-recorded  onto  diaital  taDe  in  th- 


Near  Shore  Laboratory  of  tie  Naval  Postgraduate  School.  This 
procedure  is  discussed  in  r?.ef.  3],  The  output  of  -his 
re-recording  is  approximately  ninty  minutes  of  real  time 
data,  from  either  the  sea  or  land  site,  on  a  digital  -ape 
with  a  density  of  300  bits  per  inch  (5PI)  .  The  storage 
density  proved  a  great  handicap  ir.  the  rapid  expansion  of 
the  program.  The  Church  Computer  Center  has  fifteen  tape 
drives,  but  only  or.e  of  these  is  capable  cf  handling  a  300 
BPI  tape.  The  competition  for  this  one  tape  drive  l^d  to 


processing  ti 
mane  ar.y  programming 


lost.  An  attempt  was  male  no  transfer  selected  tapes  into 
the  memory  of  the  computer,  but  it  was  unsuccessful  due  to 
miner  recording  errors  on  the  digital  tapes.  These 
recording  errors  were  overcome  by  the  read  subroutine  in 
order  to  be  able  to  work  with  the  lata.  Time  did  not  permit 
the  presuit  of  this  avenue  of  improvement  but  at  best  only  a 
few  tapes  could  be  stored  this  way  because  of  the  vast 
amount  of  raw  data. 

The  computer  program  ^as  validated  through  the  use  of  a 
test  tape.  The  test  taoa  was  recorded  with  a  sinusoidal 
input  at  known  frequency  and  voltage  level  as  measured  by  a 
sosetrum  analyser.  When  this  tape  was  processed  by  the 
landxyz  program,  the  correct  frequency  and  amplitude  were 
produced,  [Ref.  11]. 

A  noise  floor  tape  was  produced  by  replacing  the  sensing 
coils  with  120  ohm  resistors.  This  tape  was  then  evaluated 
2y  the  Lar.dtyz  program  to  produce  the  PSD  cf  the  noise  floor 
of  the  system.  figure  3.25  shows  this  noise  floor  superim¬ 
posed  or.  a  land  PSD  plot. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


.  CONCLUSIONS 


s  a  rl  i  e  r 


was  to  calibrate,  install  sr.d  operate  the  lar.i  portion 
the  Naval  Postgraduate  School’s  creomaanetic  data  ccllecti 


sys-em.  This  was  generally  successfully  e compile  had  with  rh 
exception  of  degradation  of  data  produced  by  using  the  nor 
r.oisy  La  Mesa  site  rather  than  the  nore  remote  Chew’s  Pidg 
site.  Exactly  how  much  of  an  improvement  in  the  qualify  o_ 
the  data  will  result  from  the  move  to  Chew's  Ridge  is  diffi 
cult  tc  predict.  It  could  possibly  diminish  the  SO  Hz  compo¬ 
nent  so  that  data  could  be  obtained  in  the  region  of  a  Hz 
Overall  some  70  hours  of  land  data  was  collected  during  -hi 
research,  of  which  approximately  twenty  hours  has  been  evnl 


u  at  e  a 


late. 


For  the  thesis's  second  objective,  the  production  of 
computer  program  and  procedures  to  interpret  the  data,  rh 
expansion  of  the  basic  program  to  compute  the  coherences 
Stokes  parameters,  degree  of  polarization,  elliptical  orien- 
taicn  angle,  and  the  ellioticity  would  appear  to  have  satis¬ 
fied  the  requirement.  The  major  drawback  to  this  system  i~ 
the  limited  amount  of  data  that  can  be  processed  in 
reasonable  time.  At  most  the  current  programs  are  capable  o 
evaluating  only  5-6  hours  of  land  and  sea  data  in  a  twenty- 
four  period.  This  is  primarily  due  to  ’■ape  drive  avail¬ 
ability,  but  the  requirement  of  the  man-hours  for  decoding 
re-recording  on  digital  tape  and  filing  these  tapes  at  th~ 
computer  center  cannot  be  discounted.  The  data  prccessin 
system  has  reached  the  stage  where  large  amounts  of  data  ca~ 
be  evaluated  but  it  is  not  yet  read/  to  operate  or.  a  ccnti 
r.uous  basis. 
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B.  RECOMMEND ATI  OHS 

The  first  r ecommendat ion  is  that  the  lata  transmission 
problem  be  solved  so  that  the  Ohew’s  ridge  sits  ran  be  oper¬ 
ated  on  a  continuous  basis.  This  has  two  advantages;  in  the 
improvement  in  the  quality  of  the  sampling  environment  and 
in  giving  a  contiuous  data  stream  that  can  be  recorded  at 
aopropriate  intervals  to  provide  a  Long  term  data  base  from 
which  lew  frequency  trends  could  be  observed. 


Improvement  coma.  ax=o 


processing.  Spe  cir  icax^y, 


be  mads  in 


i r e a  or 


a  new  method  of  data  handiir.a  i. 


needed.  Several  alternatives  exist,  the  first  being  the 


aeg  uisiticn 


the  projects  own  decoder  and  digital 


drive  to  previde  a  digital  tape  of  a  higher  density  than  can 
now  be  made.  This  would  make  it  possible  for  a  two  *0  eight 
fold  reduction  in  the  numoer  cf  tapes  necessary  to  hold  the 
data  and  would  also  give  faster  tarn  around  time  during 
processing  due  to  the  increased  number  of  tape  drives  avai¬ 
lable  at  the  computer  center.  Other  data  handling  schemes 


.so  feasible  if  toe  problem  of  how  -< 


aata  from 


the  digital  tapes  into  the  mass  storage  system  of  the 
computer  can  be  overcome.  This  would  make  it  possible  tc  run 
programs  in  other  than  the  ’’batch"  mode  and  avoid  the  use  of 
tape  drives  to  gain  access  to  the  data. 
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APPENDIX  A 


COiiPJ  TER  PROGRAM 

//LANDXYZ  J03  (1  573,0129),'  P03UE  SMC  1573'  , CLA  SS  =  F 
//♦MAIN  OP.G=NPG  V  Ml .  1  573P,LI  NES=(99» 

//♦FORMAT  PR, DON AME  =  PLOr .3YSVECTR,DE3  T= LOCAL 
//  EXEC  FRTXCL3  P  , REGION.  GD=  300GK 
//FORT.  SYSIN  DD  ♦ 

I  NT  EGER  *2  IN  (16) 

C  ARRAY  'IN'  IS  USED  IN  READING  DATA  FROM  TAPE 

CO  M?LEX*8  XX  ( 8  1  92)  ,  Y  Y  (3192)  ,2  (3  192) 

CGMPLEX*9  SXY  (8192)  ,  SYZ  (8  1921  ,  S3  X  (819  2) 

C  THE  COMPL  EX*  9  ARRAYS  ARE  USED  TO  ORDER  INPUT 

C  DATA  AND  INITIALLY  REPRESENTS  VOLTAGE-TIME 

C  SERIES  INFORMATION 

DIMENSION  TIME  (819  2)  ,  FREQ  (8192!  ,  WORK  (16  384)  ,  FRQ2  (91  92) 
DIMENSION  ZX  (31  92)  ,2Y  (8192)  ,ZZ(3192) 

DIMENSION  SOXY  (8192)  ,  SOYS  (8192)  ,S02X  (8192) 

DIMENSION  S1XY  (319  2)  ,  SI  YZ  (  S 1 9 2 1  ,  S1ZX  (8192) 

DIMENSION  S2XY  (8192)  ,  S2  YZ  (8192)  ,S2ZX  (8192) 

DIMENSION  S3XY  (8192)  ,  S3YZ  (8192)  ,  S3ZX  (8192) 

DIMENSION  COXY  (8 1 9  2)  ,  CO  YZ  (  8  1  9  2)  ,  COZX  (3192) 

DIMENSION  AXY  (8192)  ,  A  YZ  (9192)  ,AZX  (8192) 

DIMENSION  PXY  (8192)  ,  ?  YZ  (8  192)  ,?ZX  (819  2) 

DIMENSION  DEXY  (819  2)  ,  DEYZ  (8192)  ,DEZX  (3192) 

C  THE  »Z»  ARRAYS  REPRESENT  FREQUENCY  DOMAIN  (FFT) 

C  MAGNITUDE  DATA  AND  ARE  EVENTUALLY  CONVERTED  TO 

C  PCWEP  SPECTRAL  DENSITY  INFORMATION.  COXY  IS  COHERENCE 

C  MEASURED  X  TO  Y.  AXY  IS  THE  ANGLE  BETWEEN  THE  MAJOR 

C  AXIS  OF  THE  ELLIPSE  AND  THE  X  AXIS.  PXY  IS  THE  DEGREE 

C  OF  POLARIZATION  IN  THE  X-Y  PLANE. 

C  DEXY  IS  THE  ELLIPTI C ITY  IN  THE  X-Y  PLANE. 


I  NT  EGER  *4  IT3  (  1 2)  / 1 2 *  0/ 

REAL*1*  RT3  (28) /29*0.0/ 

REAL  ALAB  (  11)  /'  X'  ,  •  Y*  ,  •  Z‘  , 1  C’  ,  ‘  S  0  •  ,  •  S  1  •  ,  •  S  2  •  ,  •  S  3  ‘  ,  •  A  * 

,'P‘, '  E  V 

REAL* 8  TITLE  (12) 

EQOIV  ALEN CE (TITLE  (1»  ,  RT3(5) ) 

ARRAYS  ' ITB' , •RTB,f • ALAB* , AID  ‘TITLE'  ARE  USED  IN 
THE  VERTS  ATEC  PLOTTER  OUTPUT. 

DATA  XX  ,  Y  Y  /IS  3 8 4*  ( 0.  ,  0. )  / 

DATA  Z/S1 9  2* (0 .  ,0.) / 

DATA  ZX,ZY/16384*0./ 

DATA  COXY,  SOXY/16384*  0.  / 

DATA  52XY,  S3XY/163  84  *  0 .  / 

DATA  S1XY,  S17Z/16384*0./ 

DATA  COYZ,  S07Z/16334*0./ 

DATA  S2YZ,  S3  Y  Z/ 16  3  84  *  0.  / 

DATA  CO  Z7  ,  SOZX /16  3  84  *  0.  / 

DATA  S23X, S32X/16334* 0./ 

DATA  ZZ,?REQ/1o384*3. / 

DATA  AXY,  AYZ/1  6384*3.  / 

DATA  AZX,S  1ZX/ 1638 4 *3  . / 

DATA  PXY,  PYZ/16384*3.  / 

DATA  PZ X, D EXY/ 16384* 3 ./ 

DATA  DEZX, DEYZ/16384* 0./ 

PI2=6 .2831 853 

THE  FOLLOWING  SECTION  READS  THE  DECODED  PCM  DATA  OFF 
THE  DIGITAL  TAPE  AND  FOR  EACH  CHANNEL  SKIPS  EVERY 
OTHER  SAMPLE  POINT. 

I SEC=  3 

IT  L= I  SEC* 3 2 

ITL=300 

DO  55  JJ=1  , ITL 

CALL  RD  (20  ,IN  ,  200,  IRE  C, IR  R) 


55  CONTINUE 


0 

I 

I 

IS 

I 

1 


s< 


l. 


i 


IFRAME=819  2 
NR  =  16 

FNR=FLO AT  (NR) 

DO  70  L  1=  1  ,NR 

C  THE  DO  LOOP  ENDING  WITH  STATES  ENT  70  ENABLES  THE  PRO- 

C  GRAM  TO  PROCESS  A  LARGE  AMOUNT  OF  DATA  BY  REPEATING 

C  THE  PROCESS  IN  BLOCS  S .  THE  DATA  POINTS  FROM  EACH  RON 

C  THROUGH  THE  DO  LOOP  ARE  ADDED  TOGETHER  AND  EVENTUALLY 

C  AVERAGED  BY  THE  THE  NUMBER  OF  TIME  THROUGH  THE  LOOP. 

C  ’NR'  REPRESENTS  THE  NUMBER  OF  DATA  SEQUENCES  TO  BE 

C  AVERAGED 

C  1  SEQUENCE  CURRENTLY  EQUALS  8192  DATA  POINTS  FOR  EACH 

C  CHANNEL  OR  256  SECONDS  AT  A  SAMPLE  RATE  OF  32  CPS. 

DC  6  0  J J= 1  rIFR  AME 

CALL  RD  (20  ,IN,  1000  ,IR  EC,IRR) 

XX  (JJ)  =IN  (2) 

YY  (J  J)  =  IN  (3) 

3  (JJ)  =  I N  ( 4  ) 
cC  CONTINUE 

WRITS  (6,2  0  0)  I  RF.  ,IR3C 

2  00  JORMAT(10X,»IRP=»,I5,  5X,’  IREE  =  '  ,16,/) 

C  THIS  SECTION  GENERATES  THE  TIME  AND  FREQUENCY 

C  ARRAYS  AND  NORMALIZES  THE  INPUT  PCM  DATA  TO  VOLTAGE 

C  IN  PREPARATION  FOR  FAST  FOURIER  TRANSFORM  TO  THE 

C  FREQUENCY  DOMAIN. 

N=8 19  2 
FN=FLOAT(N) 

DELTAT=  1.  /  32. 

T=FN*DELT A  T 
DELTAF=  1.  0/T 
DO  20  J  =  1 , N 

TIME(J)  =D  E  LI  AT  *FLO  AT  (  J) 

FREQ  (J)  =DELTAF*FLOAT(  J) 

XX  (J)  =(  XX  (J)  -2  045. 5)  *  10. /2  ou  3.  5 
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XX  (J)  =REA  L  (XX  ( J)  ) 

YY  (J)  =  (YY  (J)  -2  045.  5|  *  10./2  04  5.  5 
YY  (J)  =REAL  (YY  (J)) 

Z(J)  =  (2  (J)  -2045.5)  *1  0  ./2D45.  5 
Z  ( J)  =  REAL  (Z(J)  ) 

'XX*  IS  THE  X-COIL  DATA,  ' YY’  IS  THE  Y-COIL  DATA, 

20  CONTINUE 

WRITE  (6,90  0)  X , Y, Z? 

01  CONTINUE 
DO  21  J  =  1  ,  N 

? HQ2  ( J)  =  ALOG  1  0  (FRSQ(J)) 

21  CONTINUE 

THE  NEXT  3  STATEMENTS  PERFORM  AN  FFT  ON  THE  INPUT 
TIME  SERIES  DATA. 

CALL  FOORT  (XX,  N, 1,-1,  0,’rfDRK) 

CALL  70URT  (YY,N,1,-1,  0,HOPK) 

CALL  FOUR!  (Z,N,  1,-1,  0  ,  W OR \  ) 

THIS  SLOCK  OF  STATEMENTS  APPLY  THE  SYSTEM  (VOLTAGE 
E-FIELD)  TRANSFER  FUNCTION  TO  THE  TRANSFORMED 
FREQUENCY  DOMAIN  DATA.  THIS  BLOCK  ENDS  AT  9. 

THE  TRANSFER  FUNCTION  CONVERTS  VOLTS  TO  NANOTESLAS 
*****WARNING***** 

THIS  TRANSFER  FUNCTION  IS  GOOD  FOR  MAGNITUDE  ONLY. 
DO  9  L=  1, N 
FFQ=FRSQ(L) 

IF  (FRQ.LE  .  25.  )  GO  IO  1 
XX  (L)  =XX(L)/28. 

GO  TO  3 

1  IF  (FRQ.  LE.  15.)  GO  TO  2 

XX(L)  =XX  (L)/(1  05.5-3.  14*FRQ) 

YY  (L)  =YY  (L)/(181. 32-7. 588*  FRQ) 

Z  (L)=Z(L)  /  (177. 26-7.  4  84*FRQ) 

GC  TO  3 

2  IF  (FRQ.  LE.  10.  )  GO  TO  3 


XX  (L)=XX(L)/(5.958*FaQ-30.97) 

YY  (L)  =YY  (L )/ (7. 166 *FRQ- 39.991 
2  (L)  =  Z  (L)  /  (6.4  9*FRQ-  3  2.35) 

GO  TO  8 

3  IF  (FRQ.  LE  .  7.  5)  GO  TO  4 

XX  (L)  =XX  (L)/(3.492*?RQ~6.  31) 

YY  (L)  =YY  (L)/(4.252*FR  Q-10.  35i 
Z  (L)  =  Z  (L)  /  (4.  044*F5Q“  7. 89) 

GO  TO  9 

4  IF  (FRQ.  LZ  .  5.)  GO  TO  5 

XX (L) =XX (L )/ ( 2 . 63 1 1  +  PRQ+O.  1466  7) 

YY  (L)  =YY  (L)/(  3.012*FRQ-1.  55) 

Z  (L)  -Z  (L)  /  (3.  184*FRQ-  1.44) 

GO  TO  8 

5  IF  (FRQ.  LE .  3.)  GO  TO  5 

XX  (L)  =XX  (L)/(2 . 63 1  1*?RQO.  1456  7  ) 

7  YY  (L)=YY  (1)/(2.702*?RQ) 

2  (L)  =  Z(L1  /  (2.  9 2*FRQ) 

GO  TO  3 

6  XX(L) =XX(L)/(2.72*FRQ) 

GO  TO  7 

8  CONTINUE 

9  CONTINUE 

THIS  BLOCK,  ENDING  A  ITH  30,  THE  MAGNITUDE  OF  THE 
FREQUENCY  DOMAIN  DATA  IS  TAKEN  AND  PLACED  IN  A  REAL 
ARRAY, 1  Z  X  ’  FOR  THE  X-COIL,'ZY’  FOR  THE  Y-COIL  AND  ZZ 
FOR  THE  Z -COIL.  THIS  DATA  IS  DIVIDED  BY  THE  NUMBER  DF 
SAMPLE  POINTS  AND  SQUARED  TD  DETERMINE  POKER. 

DO  30  11=  1  ,N 

ZX  (1 1 )  =  ZX  (II)  ♦  (CABS  (XX  (II)  )  /FN|  **2 

ZY  (II)  =ZY  (II)  ♦  (CABS  (YY  (II)  )  /FNt  **2 

ZZ  (II)  =ZZ  (II)  ♦  (CAB  S  ( Z  (II)  )  /FN)  **2 

SXY  (II)  =SXY(I1)  ♦  ((XX  (II)  /  FN)  *C0NJG  (YY  (II)  /FN)  ) 

COXY  (II)  =COXY  (II)  *CA3S(  (XX  (II)  /FN)  *CONJG  (YY  (II)  /FN)  ) 


SYZ  (II)  =S  Y  Z (1 1 }  +  ( <  YY  ( II)  /FN)  *CDNJG  (Z (I 1 ) /FN) ) 

COYZ  (II)  =CCYZ  (II)  +  CAB  S  ( (Y  Y(I1)  /?N)  *CONJG  (Z  (1 1)  /FN)  ) 

SZX  (II)  =SZX(I1)  +( (Z(I  1)  /FN)  *CON  JG  (XX  (II)  /FN)  ) 

COZX  (II)  =  COZX  (II)  4-CABS  (  (Z  (II)  /FN)  *CONJG  (XX  (II)  /FN)  ) 

30  CONTINUE 
70  CONTINUE 

THIS  BLOCK  AVERAGES  THE  DATA  POINTS  ADDED  IN  BLOCK  30 
ABOVE  BY  'NR'.  AT  THIS  POINT  THE  POWER  CURVES  ARE 
CONVERTED  INTO  POWER  SPECTRAL  DENSITIES. 

THE  COHERENCES  AND  STOKES  PARAMETERS  ARE 
COMPUTED  IN  ALL  THREE  PLANES. 

DO  33  13=1  ,N 

S2XY  (T3)=  (REAL  (SXY  (13))  /FN?)  *4. 

S3XY  (13)  =  (AIMAG  (SXY(I  3)  )  /FN?)  *4. 

S2YZ  (T3)  =  (  REAL  (SYZ  (13  )  )  /FNR)  *4. 

S3 YZ  (13 )  =  (AIMAG  (SYZ  (I  3)  )  /FNR)  *4  . 

S2EX  (13)  =  (REAL  (SZX  (1 3  )  )  /?  NR)  *4  . 

S3ZX(I3)  *  (  AIM  AG  (SZX  (I  3)  )  /  FNR)  *4  . 

COXY  (13  )=  (COXY  (13)  )  /3QRT(ZX(l3t  *ZY  (T3)  ) 

COYZ  (13)  =  (COYZ  (13)  )  /S  QRT(  ZZ(I3)  *ZY  (13)  ) 

COZX  (13)  =  (COZX  (13)  )  /3  QHT(ZX  (13)  *ZZ  (13)  ) 

ZX  (13)  *ZX  (13)  *T/FNR 
ZY  (13)  =ZY  ( 13)  *T/FN  R 
ZZ  (13)  =ZZ  ( 13)  *T/FNR 
S1XY(I3)  =  (  (ZX  ( 1 3)  -  ZY  ( 13 ) )  /T) *  2. 

S1YZ(I3)=  (  (ZY  (13) -ZZ  (13))  /T)  *2. 

S 1 Z X  (13  )  *  (  (ZZ  (13)  -ZX  (13))  /?)  *2. 

SOXY  (13)  =  (  (ZX  (13)  *ZY  (13))  /T)  *2. 

SOYZ  (13 )=  (  (ZZ  (13)  ♦ZX  (I3)|  /T)  *2. 

S0ZX(I3)  =  (  (ZX(I3)+ZZ(I3))  /T)  *2. 

33  CONTINUE 

THE  NEXT  DO  LOOP  CONVERTS  THE  3-FIELD  POWER  SPECTRAL 
DENSITY  INTO  DECIBELS  REFERENCED  TO  1  NANOTESLA  **2 
PER  HERTZ. 


ZX  (I)  =10.  *ALOG10(ZX(I  )  ) 

ZT  (I)  =10.  *  ALOG  10  (ZY  ( I )  ) 

ZZ  (I)  =10.  *  ALOG  10(ZZ{I)) 

32  CONTINUE 

THE  NEXT  BLOCK  ELIMINATES  THE  FIRST  10  DATA  POINTS  SO 
THAT  THE  FREQUENCY  IS  .04  TO  10  HZ. 

S 1 , S2  ,  AND  S3  ARE  NORMALIZED  BY  SO. 

DEGREE  Of  POLARIZATION,  ORIENTATION  ANGLE  ALPHA, 

AND  THE  ELLIPTICITY  ARE  ALSO  COMPUTED. 

DC  950  14=  1,81  32 
FREQ  (14 )  =  ?  REQ  (14+10) 

7RQ2  (14)  =F?Q2  (14+10) 

ZX  (14)  =ZX  (14+10) 

Z  Y  (14)  =Z7  (14+10) 

ZZ  (I4)=ZZ  (14+10) 

COXY  ( 14)  =CCXY  (14+1  J) 

S1XY(T4)=S  1XY  (14+10)  /SOXY  (14+10) 

S  2ZY  (14)  =  S  2XY  (14+1  0)  /  SOXY  (14+10) 

S3XY  (14 )  =  S  3XY  (14+10)  /  SOXY  (14+10) 

PXY  (14)  =SQ  RT  ( S  1  XY  ( 14 1  **2+S2XY  (14)  **2  +  S3 XY  ( 14  )  ** 2) 

DEXY  (14)  =0  .5*ARSIN  (S 3  XY  (14 )  /? X  Y  (14  )  ) 

SOXY  (14)  =  1  0.*  A  LOG  10  (SOXY  (1 4+101  ) 

COYZ  (14)  =C  OYZ  (14+10) 

S1YZ  (14)  =  S  1YS  (14+1  0)  /  SO  YZ  (14+10) 

S 2YZ  (14)  =  S  2YZ  (14+1  0)  /SOYZ  (14  +  10) 

S3YZ  (14)  =  S  3YZ  (14+10)  /  SOYZ  (14  +  10) 

P  YZ  (14)  =SQRT  (SI  YZ(I4)  **2+S2YZ  (14)  **2+S3YZ  (14  )  **2) 

DEYZ  (14)  =0.5*ARSIN  (S3  YZ  (14)  /?YZ  (14)  ) 

SOYZ  (14)  =  1  O.  +  ALOGIO  (SOYZ  (14+10)  ) 

COZX  ( 14)  =COZX  (14+10) 

SI  ZX  (14)  =S  1ZX  (14+10)  /SOZX  (14+10) 

S2ZX(I4)  =S  2ZX  (14+10)  /  SOZX  (14+10) 

S3ZX(I4)  =  S3ZX  (14+10)  /  SOZX  (14+10) 


PZX  (14)  =SQRT  (S  1ZX(I4|  **2+S2ZX  (14)  **2*S3ZX  (14)  **2) 
DEZX  (14)  =  0  .5*  A  ESIN  (3  3  ZX  (I  4)  /PZX  (1 4)  ) 

S0ZX(I4)  =  1  0.*ALO3l  0  (SOZX(I4+10)  ) 

AXY  (14)  =AT  AN  (S2XY  (14)  /S1XY(I4)  t  /2  . 

AYZ  (14)  *A  T  AN (S  2YZ ( 14 )  /S1YZ  (14)  )  /2. 

ASX  (14)  =AT  AN  (S2ZX  (I4|  /S1ZX(I4)  |  /2. 

950  CONTINUE 
C 

C  VI?  SAT EC  PLOTS 

C 

N?T3=10./0  ZLTAF  +  1  . 

C  1  NPTS  *  DETERMINES  NUMBER  OF  POINTS  NECESSARY 

C  FOR  THE  0  TO  10  HERTZ  RANGE  TO  BE  PLOTTED. 

C  FOR  THE  'ITS'  AND  »  RT3'  VALUES  REVIEW  THE  WRITE-UP 

C  FOR  THE  SUBROUTINE  PROCEDURE  'DRAW?*. 

IT3  (3)  =7 
ITB  (41=5 
IT3  (7)  =1 
ITS  (12)  =1 
RTE  (1)  =0.0 
BTE  (2)  =0.  0 
PTB  (3)  =AL  A  B  (1 ) 

READ(5,3000)  TITLE 

CALL  DRAW?  (NPTS, FRQ2,  ZX, ITS, RT3) 

RTE  (3)  =ALA  3(2) 

READ  (5,  30  00)  TITLE 

CALL  DRAWP  (NPTS, FRQ2,  ZY, ITB,  RT3) 

RTB  (3)  =AL  A  B (3) 

READ  (5, 3000)  TITLE 

CALL  DRAWP  (NPTS, FRQ2,  ZZ, I  TB,RT3) 

RTB  (3)=ALAE(4) 

READ(5,  30  0  0)  TI  TLE 

CALL  DRAW?  (NPTS,  FRQ2,  COXY  ,  IT3,  3  TB) 

RTB  (3)  =  AL  A  3  (5 ) 


77 


READ  (5,  30 0  0)  TITLE 

CALL  DRAM?  (N?TS,FRQ2,  SOXY  ,  ITS,  RTB) 
RTB  (3)  =ALA  B (6 ) 

READ  (5,3000)  TITLE 

CALL  DRAW?  (N?TS,FRQ2,  Si  XY  ,  ITB ,  R  T3) 
RTB  (3  )  =  AL  A  B  (7) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (N?TS ,  FRQ2,  S2XY  ,  ITB,  RTB) 
P.TB  (3)  =  ALA  2  (3) 

READ  (5,30  0  0)  TITLE 

CALL  DRAW?  (N?TS,FRQ2,  S3XY  ,IT3,RTB) 
RTB  (3)  =ALAE(U) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (N?T  S,  FRQ2,  COY2  ,  ITB,  R  TB) 
RTB  (3)  =ALA  B(5) 

READ  (5, 3000)  TITLE 

CALL  DRAW?  (NPT  S  ,  FR  2  ,  SO  YZ  ,  I T3  ,  ?  T  5) 
RTB  (3)  =  AL  A  5  (5 ) 

READ  (5, 30  0  0)  TITLE 

CALL  DRAW?  (N?7  S  ,  FRQ2 ,  SUE  ,  IT3,  RTB) 
F.TE  (3)  =AL  A  B  (7) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (N?TS,FRQ2,  S2  YZ  ,  IT3  ,  R  TB ) 
RTB  (3)  =ALAB(8) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS,FRQ2,  S3YZ  ,  ITB,  R  13) 
RTB  (3)  =AL  A  3  (4 ) 

READ  (5, 30  0  0)  TITLE 

CALL  DRAW?  (N?TS  ,FRQ2  ,  COZX  , ITB,  RTB) 
RTB  (3)  =  ALA  B  (5) 

READ  ( 5,  3000)  TITLE 

CALL  DRAW?  (NPTS,  FRQ2,  SOZX  ,  ITB,  3  TB) 
RTB  (3)  =ALA  B(6) 

READ  (5, 3000)  TITLE 


* 


CALL  DRAW  P  (NPTS,  FRQ2,  S1ZX,  ITB,  R  T9) 
RTE  (3)  =ALA  B  (7) 

READ  (5,  30  0  0)  TILLS 

CALL  DRAW?  (NPTS,FRQ2,  S2ZX  ,  IT3,  5  T3) 
RTB  (3)  =ALA  3(8) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS  ,FRQ2  ,  S3ZX  ,  IT3  ,  S  TB) 
RTB  (3)  =  ALA  E  (9) 

READ  (5,  30  0  3)  TITLE 

CALL  DRAW  p  (N?TS,?R02r  AX  Y,  ITB,  RTB) 
RTB  (3)  =  ALA  3(10) 

READ  (5,  3000)  TITLE 

CALL  DRAW?  (N?TS,FRQ2,  ?XY,IT3,PT3) 
RTB  (3 )  =  AL  A  H  ( 1  1) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (N?TS,FPQ2,  DEXY  ,  ITB,  STB) 
:  TB  (3)  *  ALA  2(9) 

READ  (5,30  0  0)  TITLE 

CALL  DRAW?  (N?T  E  ,?  =  52  ,  AYE  ,  IT2  ,  ?T  3) 
RTB  (3)  =  ALA  5(10) 

READ  (5,  3000)  TITLE 

CALL  Da  AW?  (N?T  S ,  FRQ2  ,  PYZ,IT3,ET3) 
RTB  (3)  =AL A  B  (1  1) 

READ  (5,3000)  TITLE 

CALL  DRAW?  (N?T5,FRQ2  ,  DEYZ  ,IT3,  RTB) 
RTB  (3)  =  ALA  3(9) 

READ  (5, 300C)  TITLE 

CALL  DRAW?  (NPTS ,  FRQ2  ,  A3X,ITB,RT3) 
RTE  (3)  =ALA  B(1  0) 

READ  (5, 30  0  0)  TITLE 

CALL  DRAWP  (NPTS  ,FRQ2  ,  ?ZX,IT3,RT3) 
RTB  (3)  =  ALAe(1  1) 

READ  (5,  30  00)  TITLE 

CALL  DRAW?  (NPTS,  FEQ2,  DEZX  ,  ITS,  RTE) 
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h: 


3300  FORMAT (6A8 ) 

STOP 
END 

//LKED.  SYSUT1  DD  S?  ACE=  ( 1  3  2  4  ,  (2  00 , 2  3  3  )  ) 
//30.SYSIN  DD  * 

LA  MESA  SITE,  16  SEP  82,  1250-1330  LOCAL 
PSD  OF  X-COIL.  AMP  IN  Dd  REF  NT**2/HZ  VS 
LA  MESA  SITE,  16  SEP  32,  1250-1330  LOCAL 
PSD  OF  Y-COTL.  AM?  IN  D3  ?.3F  NT**2/HZ  VS  HZ 
LA  ME  S A  SITE,  16  SEP  82,  12  50-1  330  LOCAL 
PSD  OF  Z-COIL.  AMP  IN  DB  ?.SF  NT**2/H2  VS  HZ 
LA  MESA  SITE,  16  SEP  32,  12  50-1  330  LOCAL 
CCHESENCS  OF  X  5  Y  COILS  VS  FREQ  (HZ) 

LA  MESA  SITE,  16  SEP  82,  12  50-1  330  LOCAL 
SO  (  10  LOG)  VS  FREQ(HZ)  X-Y  PLANS 

LA  MESA  SITE,  16  SEP  82,  12  50-1  330  LOCAL 
SI/SO  7 S  FREQ  (HZ)  X-Y  PLANE 

LA  MESA  SITE,  16  SEP  32,  1250-1330  LOCAL 

52/SO  VS  FREQ  (HZ)  X-Y  PLANE 

LA  MESA  SITE,  16  SEP  82,  1  2  50-1  330  LOCAL 
S3/S0  VS  FREQ  (HZ)  X-Y  PLANE 

LA  MESA  SITE,  16  SEP  82,  1250-1330  LOCAL 
COHERENCE  OF  Y  6  Z  COILS  VS  FP.SQ(HZ) 

LA  MESA  SITE,  16  SEP  82,  1  2  50-1  3  30  LOCAL 
SO  (  10  LOG)  VS  FREQ(HZ)  Y-Z  PLANS 

LA  MESA  SITS,  16  SEP  82,  12  50-1  330  LOCAL 
Si/30  VS  FRSQ(HZ)  Y-Z  PLANE 

LA  MESA  SITE,  16  SEP  32,  12  50-1330  LOCAL 
S2/S0  VS  FREQ  (HZ)  Y-Z  PLANE 

LA  MESA  SITE,  16  SEP  82,  12  50-1  330  LOCAL 
S3/S0  7 S  FREQ  (HZ)  Y-Z  PLANE 

LA  MESA  SITS,  16  SEP  82,  12  50-1  330  LOCAL 
COHERENCE  OF  Z  5  X  COILS  VS  FREQ  (HZ) 

LA  MESA  SITE,  16  3Z?  82,  1250-1330  LOCAL 
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SO  ( 1  0  LOG )  VS  FESQ(HZ) 

z-x 

PLANE 

LA  MESA 

SITE,  16  SEP 

82, 

12  50-1  330 

LOCAL 

SI/SO  vs 

FREQ  (HZ) 

Z-X 

PLANE 

LA  MESA 

SITS,  16  SEP 

82, 

1250-1  330 

LOCAL 

S 2/SO  VS 

FREQ  (HZ) 

Z-X 

PLANS 

LA  MESA 

SITE,  16  SEP 

82, 

1  2  50-1  330 

LOCAL 

S3/S0  VS 

FREQ  (HZ) 

Z-X 

PLANE 

LA  MESA 

SITE,  16  SEP 

82, 

12  50-1  330 

LOCAL 

A.i’i  7  S  F 

ESQ  (HZ) 

LA  MESA 

SITE,  16  SEP 

82, 

1  2  50-1  3  30 

LOCAL 

DEG~.EE  CF  POLARIZATION  VS 

F  REQ(K  Z) 

X-Y 

LA  MESA 

SITE,  16  SEP 

82, 

1  2  50-1  330 

LOCAL 

7f ” Twr^T 

Li-Lj  lLl 

TY  VS  FREQ  (HZ) 

X-Y 

PLANS 

LA  MESA 

SIT E,  16  SEP 

82, 

1  2  50-1  330 

LOCAL 

AYE  V  Sr 

EEQ  (HZ) 

LA  MESA 

SITE,  16  SEP 

82, 

12  50-1  330 

LOCAL 

J  EG  P.  E  E  C 

F  POLARIZATION  VS 

FREQ  (HZ) 

Y-Z 

LA  MESA 

SITE,  16  SEP 

82, 

12  50-1  3  30 

LOCAL 

ELL ITICI 

TY  VS  FREQ  (HZ 

;) 

Y-Z 

PLANE 

LA  MESA 

SITE,  16  SEP 

8  2, 

1  2  50-1  3  30 

LOCAL 

A  EX  VS  E 

ESQ  (HZ) 

LA  MESA 

SITS,  16  SEP 

82, 

12  50-1  330 

LOCAL 

DEGREE  CF  POLARIZATION  VS 

F  REQ(HZ) 

Z-X 

LA  MESA 

SITE,  16  SEP 

82, 

1  2  50-1  330 

LOCAL 

ELL ITICI 

TY  VS  FREQ  (HZ 

:) 

Z-X 

PLANE 

/* 

//SC. FT20F001  DD  UN IT=3U 03 - 4 , VO L=SER=S K DT52 , DIS P= (OL D , KEEP)  , 
//  LA  BSL  =  (1  ,NL,  ,  IN)  , 

//  DC  E=  (RECFM=F3,LRECL=32,  3  L  KSIZ  S=  5  1  2  ,  D  SN  =  2) 

/* 

// 
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APPENDIX  B 


SJJISIS  ri d j s  sirs  PREPARATION 

Chew's  Ridqe  is  the  sire  of  previous  land-side 
geomagnetic  data  collection.  It  is  located  approximately 
twenty-three  air  miles  southeast  of  the  Naval  Postgraduate 
School  or.  a  hilltop  3,303  fe=t  Verve  =  l  =  v  =  l.  It  was 
chosen  for  its  remoteness  from  -he  local  power  grit  (3-4 
miles)  and  its  virtual  "line  of  sight"  RF  path  to  the 
school.  A  detailed  description  of  the  physical  layout  an! 
operation  of  the  site  car.  be  found  in  [Ref.  5]. 

The  original  intention  of  this  research  was  to  collect 
data  from  this  remote  site  to  provide  the  simultaneous  data 
for  the  sea  trials.  Preparation  of  the  site  and  the 
necessary  equipment  continued  througa  most  of  the  project. 
The  da-a  transmission  from  this  site  to  the  Naval 
Postgraduate  School  receiving  location  proved  impossible 
with  the  transmitter/receivers  available. 

Previous  data  collected  from  Chew's  Ridge  was  collected 
using  a  Motorola  system  that  suffered  from  overheating 
failures  when  required  to  operate  in  a  continuous  mode. 
Intially  another  attempt  was  made  to  repair  the  Motorola 
unit  but  after  one  last  failure  this  proved  economically 
impossible  and  it  was  felt  that  this  unit  did  not  have  the 
bandwidth  to  accomadate  the  addition  of  the  PCM  to  the 
system.  Cor.segu  ently  a  RSPCO  Model  310  transmit  t  er- r  eceiver 
pair  was  ordered.  These  on  arrived  the  first  of  September 
and  initial  testing  showed  that,  although  they  were  capable 
of  operating  for  long  periods  of  time,  they  did  not  have  the 
required  29  KHz  bandwidth  to  accomadate  the  PCM  data. 
Efforts  are  currently  underway  to  expand  the  bandwidth  to 
meat  the  requirements. 
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To  properly  orient  the  sensing  coils  at  this  site  it  is 
necessary  to  have  an  accurate  known  direction.  In  order  to 
accomplish  this,  inquiries  were  made  to  the  National  Forest. 
Service  to  determine  if  they  have  survey  data  for  the  area. 
Whan  it  proved  nonexistent,  a  survey  by  resection  was 
conducted.  A  theodolite  was  used  to  measure  relative  angles 
to  several  prominent  terrain  features  in  the  area  from  the 
Chew’s  Ridge  site.  These  measurements  were  made  three  times 

a  ~  h  *  V  pi  r  m  ^  ^  ^  3  —  *“  i  I  ^  '  73  a  3 "  a  ^  c.  ~  a 

drawn  or.  a  piece  of  overlay  taper  and  oriented  on  a  1  :  50,000 
mao  of  the  area.  This  allowed  for  the  best  fit  of  the  data 
and  prevented  blunders  in  the  recognition  of  the  terrain 
features.  The  most  prominent  terrain  feature  (Uncle  Sam 
Mountain)  was  then  chosen  for  the  end  of  the  orienting  line. 
The  Grid  Azimuth  tc  Uncle  Sam  Mtn.  was  then  computed  from 
the  coordinates  of  the  two  locations.  Finally  the  True 
Azimuth  was  calculated  to  be  25  7. '4  degrees.  Es-ima*ed 
accuracy  of  this  measurement  is  or. e  uegree. 
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